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Abstract 

We make a QCD-inspired parameterization of all accelerator data on forward proton-proton and 
antiproton-proton scattering amplitudes. Using vector dominance and the additive quark model, we 
show that the same parameters also fit yp and 77 interactions. Using the high energy predictions of our 
model, along with Glauber theory, we calculate proton-air cross sections at energies near y/s « 30 TeV. 
The comparison of p-air cosmic ray measurements with our QCD model predictions provide a strong 
constraint on the inclusive particle production cross section. 
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1 Introduction 



This communication is divided into three sections. 

First, we show that the data on the total cross section, the slope parameter B of the elastic differential 
cross section, and the ratio of the real to imaginary part of the forward scattering amplitude p for pp and 
pp interactions can be nicely described by a model where high energy cross sections grow with energy as a 
consequence of the increasing number of soft partons populating the colliding particles Q , ■ The differential 
cross sections for the Tevatron and LHC are predicted. 

Next, we verify the model by showing that the known experimental data on jp and 77 interactions can 
be derived from our pp and pp forward scattering amplitudes using vector meson dominance (VMD) and the 
additive quark model[||. 

Finally, we use the high energy predictions of our QCD-inspired parameterization of accelerator data on 
forward proton-proton and antiproton-proton scattering amplitudes, along with Glauber theory, to predict 
proton-air cross sections at energies near y/s ~ 30 TeV[[|. 

All cross sections will be computed in an eikonal formalism guaranteeing unitarity throughout: 

a tot (s) = 2 J {l-e-M^cosIx^s)]} d 2 b. (1) 

Here, x is the complex eikonal (x = Xj? +*Xj)i an d b is the impact parameter. The even eikonal profile function 
^even rece j ves contributions from quark-quark, quark-gluon and gluon-gluon interactions, and therefore 

X even (s,b) = x qq (s,b)+x qa (s,b)+x gg (s,b) 

= i [a qq (s)W(b; fj, qq ) + a qg (s)W(b; yj\x qq \x gg ) + a gg (s)W(b; fi gg )] , (2) 

where cr.y are the cross sections of the colliding partons, and W(b; fi) their overlap function in impact 
parameter space, parameterized as the Fourier transform of a dipole form factor. The impact parameter 
space distribution function 

W(b;v) = ^-(pb) 3 K 3 (pb) (3) 
967T 

is normalized so that J W(b ; fi)d 2 b = 1. As a consequence of both factorization and the normalization chosen 
for the W(b ; /it), it should be noted that 

X even (s, b) d 2 b = 1 [a gg (s) + a qg (s) + a qq (s)\ , (4) 

so that <rfo t era (s) = 2Im{i [c gg (s) + <J qg (s) + a qq (s)]}, for small x- This formalism is identical to the one used 
in "mini-jet" models B, as well as in simulation programs for minimum-bias hadronic interactions such as 
PYTHIA and SIBYLL||. 

In this model hadrons asymptotically evolve into black disks of partons. The rising cross section, asymp- 
totically associated with gluon-gluon interactions, is simply parameterized by a normalization, an energy 
scale, and two parameters: /i gg which describes the "area" occupied by gluons in the colliding hadrons, and 
J(= 1 + e). Here, J is defined via the gluonic structure function of the proton, which is assumed to behave 
as l/x J for small x. It therefore controls the soft gluon content of the proton. The introduction of the 
quark-quark and quark-gluon terms allows us to adequately parameterize the data at all energies, since the 
"size" of quarks and gluons in the proton can be different. In the present context, this model represents a 
convenient parameterization of the pp and pp forward scattering amplitude. 

The photoproduction cross sections are calculated from this parameterization assuming vector meson 
dominance and the additive quark model. For the probability that the photon interacts as a hadron (Phad), 
we use the value Phad — 1/240 which can be derived from vector meson dominance. Our results show that 
its value is indeed independent of energy. It is, however, uncertain by 20% because it depends on whether 
we relate photoproduction to 7r-nucleon or nuclcon-nucleon data (In other words, irN and NN total cross 
sections only satisfy the additive quark model to this accuracy). Subsequently, following reference ||, we 
obtain jp cross sections from the assumption that, in the spirit of VMD, the photon is a 2 quark state in 
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contrast with the proton which is a 3 quark state. The jp total cross section is obtained from the even 
eikonal for pp and pp by the substitutions — ► | cry and fii — ► /i^ . 

We will thus produce a parameter-free description of the total photoproduction cross section, the phase 
of the forward scattering amplitude and the forward slope for 7p — > Vp, where V = p,uj, <p. Interestingly, our 
results on the phase of Vp — > Vp are in complete agreement with the values derived from Compton scattering 
results (7 + p — > 7 + p) using dispersion relations. We also calculate the total elastic and differential cross 
sections for "/p — > Vp. This wealth of data is accommodated without discrepancy. 

The 77 cross sections are derived following the same procedure. We now substitute cry — > | cry and 
Pi — > I Hi into the nucleon-nucleon even eikonal, and predict the total cross section and the differential cross 
sections for all reactions 77 — > ViVj at a variety of energies, where V — p, w, <f>. 

The high energy 77 total cross section have been measured by two experiments at LEP. While these 
measurements yield new information on its high energy behavior at center-of-mass energies in excess of 
\fs =15 GeV, they may represent the last opportunity to measure the 77 cross section, and the two data 
sets appear to disagree. However, it has been argued that the original data are consistent within the errors 
and that the observed disagreements are due to two different Monte Carlo's used to extract the quoted 
values. We here point out that our analysis nicely accommodates the L3 result J8|. Our model approximately 
satisfies the factorization theorem, o vv ja lp = cr 7 j,/cr 77 , because of its small eikonal. The OPAL data do 
not satisfy it. In fact, no model incorporating the additive quark model and factorization can accommodate 
the OPAL data. VMD and factorization are sufficient to prevent one from adjusting Phad, or any other 
parameters, to change this conclusion. 



2 High energy proton-proton and proton-antiproton scattering 

In this section we discuss our QCD-inspired parameterization of the forward amplitudes. To determine its 
parameters, wc fit all high energy forward pp and pp scattering data above 15 GeV, for the total cross section 
(ctot)) the ratio of the real to the imaginary part of the forward scattering amplitude (p), and the logarithmic 
slope of the differential elastic scattering cross section in the forward direction (B). Then, we compare the 
experimental data for the elastic scattering cross section and for the differential elastic scattering with our 
results. Finally, a prediction is made for the differential elastic scattering at the LHC. 

Our QCD-inspired parameterization satisfies crossing symmetry, i.e., it is either even or odd under the 
transformation E — ► — E, where E is the laboratory energy. This allows us to simultaneously describe 
pp and pp scattering. It also satisfies analyticity, and unitarity because of the eikonal formalism. Since 
the total cross section asymptotically rises as log 2 s, our QCD-inspired parameterization complies with the 
Froissart bound. The eikonal formalism for calculating a tot-, P and B, along with details on the analyticity, 
the Froissart bound, and the QCD-inspired eikonal are given in ref. In all 11 parameters are used. The 
low energy region, where the differences between pp and pp scattering are substantial, largely determines 
the 7 parameters necessary to fit the odd eikonal and the quark-quark and quark-gluon contribution to the 
even eikonal. Thus, they largely decouple from the high energy behavior. Hence, for y/s > 25 GeV, where 
the difference between pp and pp scattering becomes small, only 4 parameters describe all data. 

We fit all the highest energy cross section data (E710 |J, CDF [|l(J and the unpublished Tevatron value 
|pT| ), which anchor the upper end of our cross section curves. The results of the fit are shown in Fig. [l]. 
Data for p values and B are confronted with our model in Figs. I and |. 

It can be seen from those figures that we obtain a satisfactory description of all 3 quantities, for both 
pp and pp scattering. The x 2 of the fit is reasonably good (considering the large spread in some of the 
experimental data, as well as the discrepancies in the highest energy cross sections), giving a x 2 /d-f- — 1-66, 
for 75 degrees of freedom. The model splits the difference between the measurements of the total cross 
section at y^i = 1800 GeV (see Fig. (IK From Fig. g, we note that the fit to p is anchored at \fs = 550 GeV 
by the very accurate measurement |1^| of UA4/2 and passes through the E710 point |l3|j . The statistical 
uncertainty of the fitted parameters is such that at 25 GeV the cross section predictions are statistically 
uncertain to ss 1.3%, at 500 GeV are uncertain to sa 1.6%, and at 2000 GeV are uncertain to ~ 2.5%. 

In Fig. H we show the prediction for the elastic cross section along with the data for both pp and pp. The 
agreement is excellent. We note that a e i a stic is rising more sharply with energy than the total cross section 
a tot . Comparing Fig. |l| with Fig. ||, we see that the ratio of the elastic to total cross section is rising with 
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Figure 1: The total cross section for pp and pp scattering. The solid line and squares are for pp and the 
dotted line and circles are for pp. 
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Figure 2: The ratio of the real to imaginary part of the forward scattering amplitude for pp and pp scattering. 
The solid line and squares are for pp and the dotted line and circles are for pp. 



- 3 - 



20 




Vs, in GeV 



Figure 3: The nuclear slope parameter for elastic pp and pp scattering. The solid line and squares are for pp 
and the dotted line and circles are for pp. 
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Figure 4: Elastic scattering cross sections for pp and pp scattering. The solid line and squares are for pp and 
the dotted line and circles are for pp. 
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energy. The ratio is, of course, bounded by the value for the black disk |T4], ffq], i.e., 0.5, as the energy goes 
to infinity. 

Having fixed all parameters specifying our eikonal, we calculate da/dt, for various values of y/s. The 
differential cross section at the Tevatron (yfs = f800 GeV) is shown in Fig. || along with E710 |T(| data. 
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Figure 5: The elastic differential scattering cross section for the reaction pp — > pp at y/s = 1 800 GeV. The 
data points are from E710. 

The agreement over 4 decades is striking. 

Our prediction for the differential cross section at y/s — 14 TeV, the energy of the LHC, is plotted in 
Fig. H In particular, at small \t\, we predict that the curvature parameter C (da/dt cx e Bt+ct for small t; 
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Figure 6: The elastic differential scattering cross section for the reaction pp — > pp at LHC. 

see ref. [Q for details) is negative. For energies much lower than 1800 GeV, the observed curvature has 
been measured as positive. For 1800 GeV, we see from Fig. || that the curvature parameter C is compatible 
with being zero. Block and Calm [|l4[ [l|| have pointed out that the curvature is predicted to go through 
zero near the Tevatron energy and that it should become negative thereafter. Asymptotically the proton 
approaches a black disk. Its curvature is always negative 0, ^B), C = — i? 4 /192, where R is the radius of 
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the disk. Thus, the curvature has to pass through zero as the energy increases. 'Asymptopia' is the energy 
region (energies much larger than the Tevatron) where the scattering approaches that of a sharp disk. 

With the parameters we obtained from our fit, the total cross section at the LHC (14 TeV) is predicted 
to be a tot — 108.0 ± 3.4 mb, where the error is due to the statistical errors of the fitting parameters. 



3 Photon-proton reactions 

We assume that the photon behaves like a two quark system when it interacts strongly. We therefore obtain 
7P scattering amplitudes by performing the substitutions <7y — > |<7y and pi — > yf Mi m the even eikonal 
for nucleon-nucleon scattering, so that 



= i [l<J qq {s)W (b; yfiiigg) + lcr qg {s)W ( 



* gg (s)w(b 



2^99 



Using vector dominance, the photon-proton total cross section is then written as 



ad I jl 



C0S [Xr 



(M)]} 



d 2 b, 



(5) 



(6) 



where Phad is the probability that a photon interacts as a hadron. We use the value Phad — 1/240. This 
value is found by normalizing the total 7p cross section to the low energy data, and is very close to that 
derived from vector dominance, 1/249. Using fp/An = 2.2, /^/4tt = 23.6 and /|/47r = 18.4, we find 
£ v (47ra//2) = 1/249, where V = p,w,(f> (see Table XXXV, pag. 393 of Ref . [0). 

With all eikonal parameters fixed by the nucleon-nucleon data, we can now calculate oJ^ t ( s ) ■ The result is 
shown in Fig. (7| It reproduces the rising cross section for jp, using the parameters fixed by nucleon-nucleon 
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Figure 7: The total cross section for scattering. 

scattering. This prediction only uses the 9 parameters of the even eikonal, of which but 4 are important in 
the upper energy region. The accuracy of our predictions are ~ 1.5%, from the statistical uncertainty in our 
eikonal parameters. 

We next consider the 'elastic' scatterings 

7 + p -> Pvirtual +P -> P + P , 
7 + p -> Uyirtual + P -> W + P , 

7 + p — > ^virtual + P — > 4> + P ■ (J) 

Here the photon virtually transforms into a vector meson which elastically scatters off of the proton. The 
strengths of these reactions is O(a) times a strong interaction cross section. The true elastic cross section is 
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given by Compton scattering on the proton, 7 + p — > 7 + p, which we can visualize as 



7 + p -> Pvirtual + P— > p+p— >7+P> 
7 + p -> Virtual +p^UJ+p^~f+p, 

J+p -> (^virtual +P -> + P -> 7 • (8) 

It is clearly C(a 2 ) times a strong interaction cross section, and hence is much smaller than 'elastic' scattering 
of Eq. (|7|). Thus, we justify the use of Eq. (^J) to calculate the total cross section, since only reactions with 
a photon in the final state are neglected. 

We evaluate p and the slope B for the 'elastic' scattering expressed in Eq. (fj]), with p and B being the 
same for all 3 reactions. 

The dependence of p with the energy is shown in Fig. ||. Damashek and Gilman []l8| have calculated 
the p value for Compton scattering on the proton using dispersion relations, i.e., the true elastic scattering 
reaction for photon-proton scattering. We compare this calculation, the dotted line in Fig. 0, with our 
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Figure 8: Ratio of the real to imaginary part of the forward scattering amplitude for the 'elastic' reactions 
/+p — > Vi +p, where Vi is p , uj° or <jp . The dotted curve is for Compton scattering from dispersion relations 
It has been slightly displaced from the solid curve for clarity in viewing. 



prediction of p (the solid line). The agreement is so close that we had to move the two curves apart so that 
they may be viewed more clearly. 

In Fig. ^ we show our results for the slope B as a function of the energy. The available experimental 
data for 'elastic' pp and uip final states are also plotted. Again, the agreement of theory and experiment is 
very good. 

To calculate the elastic cross sections cr^ stic and differential cross sections da Vp /dt as a function of 
energy, we use 



Vp t \ 
® elastic V / 



had 



1 



d 2 L 



(9) 



where PuJj is the appropriate probability for a photon to turn into V, with V = p, lo or cf>. The differential 
scattering cross section is given by 



da y P 
dt 



(*.*) 



pVp 
had 

47T 



J (g6)(l-e IX " P(M) )d 2 -b 



(10) 



where t — —q 2 . 

Since we normalize our data to the cross section found with x 1P i an d n °t to (cr^ ot + cr^ ot )/2, we must 
multiply all fy/in by 1.65. Hence, our effective couplings are j 2 pc sl^ = 3-6, /„ e ff/^ 7r = 38.9, and /| e g/47r = 
30.4. 
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Nuclear Slope 
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■ 7+p— > co+p, Zeus 

A Y+p-> P +p, low energy 
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Figure 9: Nuclear slope parameter for the 'elastic' reaction 7 + p — > V, L + p, where Vi is p°, uj° or <f>. For 
the reaction 7 + p — > p° + p, the inverted triangles are the Zeus data, the circles are the HI data, and the 
triangles are the low energy data. For the reaction 7 + p — > w° + p, the squares are the Zeus data. 
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shown in Figs. [l(] and [ll], respectively. 
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Figure 10: The 'elastic' photoproduction cross section, for the reaction ~f+p — > p +p. The squares are Zeus 
data, the circles are HI data, and the inverted triangles the low energy data. 

The differential cross section, da/dt, for the 'elastic' reactions 7+p^ p° + p, 7 + p — > w° + p and 

, and respectively. The agreement, in absolute normalization 
and shape, of our results for all three light vector mesons with the experimental data for all available energies 
reinforces our confidence in the model. 



7 + p — > <j) + P are plotted in Figs. [12, 13 



4 Photon-Photon Interactions 

In this section, we consider 77 interactions. As it was done for jp interactions, we will start from the eikonal 
X yp {s 1 b) and multiply every cross section by 2/3 and multiply each p by y3/2. Therefore, 

X ^(s, b) = i [la qq (s)W (b; + ^a qa {s)W (b; + t^(^)^ (b; \p gg )] . (11) 
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Figure 11: The 'elastic' photoproduction cross section for the reaction 7 + p — > u>° + p. The circles are Zeus 
data, and the squares are the low energy data. 
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Figure 12: The differential cross section for the 'elastic' reaction 7 + p — > p° + p. The solid curve and the 
circles (Ballam et al. data) are at \/s— 4.3 GeV, the dashed curve and triangles (HI data) are at \fs= 55 
GeV, and the dotted curve and diamonds are at y/s= 73 GeV (Zeus data). 
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Figure 13: The differential cross section for the 'elastic' reaction 7 + p —> ui° + p at -^5=80 GeV. The circles 
arc the Zeus data. 
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Figure 14: The differential cross section for the 'elastic' reaction 7 + p ■ 
are the Zeus data. 



-p at -^/s^TO GeV. The circles 



- 10 - 



Using vector dominance we obtain, 

<T%(s)=2I# ad J {l-e^ 7 M cos[f/(M)] } d % (12) 
where Phad — 1/240 is the probability that a photon will interact as a hadron. In Fig. Il5] we plot our results 
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Figure 15: The total cross section for 77 scattering. The data sources are indicated in the legend. 

for a2g t {s) as a function of the energy, and compare it to the various sets of experimental data. We note that 
our prediction fits the L3 data, but doesn't fit the OPAL results. 

5 Proton-air cross sections 

Cosmic ray experiments measure the penetration in the atmosphere of particles with energies in excess of 
those accelerated by existing machines — interestingly, their energy range covers the energy of the Large 
Hadron Collider (LHC) and extends beyond it. However, extracting proton-proton cross sections from 
cosmic ray observations is far from straightforward pj| . By a variety of experimental techniques, cosmic 
ray experiments map the atmospheric depth at which cosmic ray initiated showers develop. The measured 
shower attenuation length (A m ) is not only sensitive to the interaction length of the protons in the atmosphere 
(Ap_air), with 

A m = fcAp_ a i r = k r— j , (13) 



but also depends on the rate at which the energy of the primary proton is dissipated into electromagnetic 
shower energy observed in the experiment. The latter effect is parameterized in Eq. ( [l3"| ) by the parameter 
k; m p is the proton mass and cr™*^ the inelastic proton-air cross section. The value of k depends on the 
inclusive particle production cross section in nucleon and meson interactions on the light nuclear target of 
the atmosphere and its energy dependence. We here ignored the fact that particles in the cosmic ray "beam" 
may be nuclei, not just protons. Experiments allow for this by omitting from their analysis showers which 
dissipate their energy high in the atmosphere, a signature that the initial energy is distributed over the 
constituents of a nucleus. 

The extraction of the pp cross section from the cosmic ray data is a two step process. First, one calculates 
the p-air total cross section from the measured inelastic cross section 

,1 _ a-el 



^inel 

°p-air = 



p — air 



p— air 



(14) 



Next, the Glauber method (20j is used to transform the measured value of off?fL r into a proton-proton total 
cross section a pp ; all the necessary steps are calculable in the theory. In Eq. (JIJ) the measured cross section 
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for particle production is supplemented with the elastic and quasi-elastic cross section, as calculated by the 
Glauber theory, to obtain the total cross section cr p _ a i r . The subsequent relation between o^fjjj and a vv 
involves the slope of the forward scattering amplitude for elastic pp scattering, 



B 



d_ 

dt 



dt 



(15) 



and is shown in Fig.|l6|, which plots B against cr pp , for 5 curves of different values of o~ 1 p n °l li] .. This summarizes 
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Figure 16: B dependence on the pp total cross section a pp . The five curves are lines of constant cr p ™\ ir , of 
440, 490, 540, 590 and 640 mb — the central value is the published Fly's Eye value, and the others are ±1<7 
and ±2cr. The dashed curve is a plot of our QCD-inspired fit of B against a pp . The dot is our value for 
yfs = 30 TeV, the Fly's Eye energy. 

the reduction procedure from cr™^ to a pp [|l9|. Also plotted in Fig. |l6|is a curve of B vs. a pp which will be 
discussed later. 

A significant drawback of the method is that one needs a model of proton-air interactions to complete 
the loop between the measured attenuation length A m and the cross section (T^ain *- e -' the va hie of k in Eq. 
(|l3|). We minimize the impact of theory by using our QCD-inspired parameterization of the forward proton- 
proton and proton-antiproton scattering amplitudes which is analytic, unitary and simultaneously fits all 
data of (Ttot, B and p. Using vector meson dominance and the additive quark models, we have shown that it 
accommodates a wealth of data on photon-proton and photon-photon interactions without the introduction 
of new parameters. Because the model is both unitary and analytic, it has high energy predictions that are 
essentially theory-independent. In particular, it also simultaneously fits o~ pp and B, forcing a relationship 
between the two. Specifically, the B vs. a vp prediction of the model is shown as the dashed curve in Fig.|l6|. 
The dot corresponds to our prediction of a pp and B at y/s = 30 TeV. It is seen to be slightly below the curve 
for 490 mb, the lower limit of the Fly's Eye measurement, which was made at yfs « 30 TeV. 

In Fig. [l7| we have plotted the values of a pp vs. gj™^ that are deduced from the intersections of the 
B-o- pp curve with the o^!^ curves of Fig. [if]. 

Figure [lj] allows the conversion of the measured Cp"^^ to (j pp 

Our prediction for the total cross section a pp as a function of energy is confronted with all of the accelerator 
and cosmic ray measurements [2l| E2[ in Fig.[lJ. For inclusion in Fig.[lJ, we have calculated the cosmic 
ray values of a pp from the published experimental values of cp 1 ^^, using the results of Fig.^j. We note the 
predicted curve is systematically lower than the cosmic ray points, roughly about the level of one standard 
deviation. 

It is at this point important to recall Eq. ( |l3| ) and consider the fact that the extraction of o^^t from 
the measurement of A m requires a determination of the parameter k. The measured depth V max at which a 
shower reaches maximum development in the atmosphere, which is the basis of the cross section measurement 
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Figure 17: A plot of the predicted total pp cross section a pp , in mb vs. the measured p-air cross section, 
mb. 



cr mcl • in 

u v— air) 111 




1e+2 1e+3 
\'s, in GeV 



1 e+4 1 e+5 



Figure 18: A plot of the QCD-inspired fit of the total nucleon-nucleon cross section a pp , in mb vs. y/s, 
in Gev. The cosmic ray data that are shown have been converted from <Tp^; r to a vp using the results of 
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in Rcf . pl| , is a combined measure of the depth of the first interaction, which is determined by the inelastic 
cross section, and of the subsequent shower development, which has to be corrected for. The position of 
Imai also directly affects the rate of shower attenuation with atmospheric depth which is the alternative 
procedure for extracting 

The model dependent rate of shower development and its fluctuations are the origin of the deviation of 
k from unity in Eq. (|l3|). Its values range from 1.5 for a model where the inclusive cross section exhibits 
Feynman scaling, to 1.1 for models with large scaling violationsp9|. The comparison between data and 
experiment in Fig. |lS| is further confused by the fact that the AGASAp2| and Fly's EyepT| experiments 
used different values of k in the analysis of their data, i.e., AGASA used k = 1.5 and Fly's Eye used k — 1.6. 

We therefore decided to match the data to our prediction and extracted a common value for k = 1.33 ± 
0.04. This neglects the possibility that k may show a weak energy dependence over the range measured. 
In Fig. [l9] we have replotted the high energy cosmic ray data for our prediction of crj^ lr vs. y/s, with the 
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Figure 19: A \ 2 fit of the measured AGASA and Fly's Eye data for u™f l ail , in mb, as a function of the 
energy -y/s, in GeV. The result of the fit for the parameter k in Eq. (|l^) is k = 1.33 ± 0.04. 

common value of 1.33 obtained from a \ 2 fit- Clearly, we have an excellent fit, with good agreement between 
AGASA and Fly's Eye. The analysis gives % 2 = 1.75 for 6 degrees of freedom (the low x 2 is probably due to 
overestimates of experimental errors). This result for k is interesting — it is close to the value of 1.2 obtained 
using the SIBYLL simulation [^4] for inclusive particle production. This represents a consistency check in 
the sense that our model for forward scattering amplitudes and SIBYLL share the same underlying physics. 
The increase of the total cross section with energy to a black disk of soft partons is the shadow of increased 
particle production which is modeled by the production of (mini)-jets in QCD. The difference between the 
k values of 1.20 and 1.33 could be understood because the experimental measurement integrates showers in 
a relatively wide energy range, which tends to increase the value of k. 

In the near term, we look forward to the possibility of repeating this analysis with the higher statistics 



of the HiRes 25 cosmic ray experiment that is currently in progress and the Auger |2(| Observatory. 

In conclusion, we have successfully united the high energy cross section results (-^s « 30 TeV) of the 
cosmic ray measurements with the accelerator cross section measurements, under a common rubric, the 
QCD-inspired analysis. 
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